High itlt pion and proton correlations in central PbPb collisions at ^/snn = 17GeV 
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For central PbPb collisions at \^S„„ = 17.3 GeV we have made the first two-dimensional measure- 
ment of the pp correlation function. These data extend the range of previous studies of HBT radii 
by a factor of two in rriT ■ They are consistent with a hydrodynamic interpretation and microscopic 
models that include hadronic rescattering and transverse expansion. We also report new data on 
pion correlations. The two particle correlations of negative pions at hit = 0.92 GeV imply source 
radii that are smaller than typical hydrodynamic fits and transport model simulations. It is possible 
that these fast pions may have left the source before the build up of hydrodynamic flow. 

PACS numbers: 25.75.Gz, 25.75.Ld 



Experimental studies of high-energy nuclear collisions 
at the AGS, SPS and RHIC accelerators (with \fS nn = 
2 — »200 GeV) have revealed many interesting features of 
hot and dense nuclear matter, and some characteristic 
signatures of a quark-gluon-plasma phase have been re- 
ported pj. If such a high density state were formed in 
the initial stages of the reaction the high initial presure 
would result in a significant transverse flow of the hadrons 
that originate from the participant zone. Depending on 
the time scale and the amount of rescattering the system 
might be expected to be in local thermodynamical equi- 
librium. This conjecture is supported by the linear in- 
crease of the single particle inverse slopes with mass H0|. 

The Hanbury-Brown Twiss, (HBT), technique uses the 
interference of particle wavefunctions to infer the an- 
gular size of stars or the length scales of nuclear sys- 
tems from the two particle correlation function Q . The 
"radii" measured by HBT can be thought of as "lengths 
of homogeneity" of the source which depend upon ve- 
locity and/or temperature gradients HH 0- Heavy 
Ion HBT measurements over a wide range of energies 
11 H El El El E l El El El show a drop of radii with 
hit = \/ Vt "h m2 consistent with a hydrodynamical in- 
terpretation [H El El ill 



Faster particles are more likely to come from initial col- 
lisions of the incoming baryons, so at some momentum we 
might expect the hydrodynamic approach to break down. 
This may have already been seen at RHIC ( v / sjvjv = 
200GeV) where the strength of elliptic flow falls below 
the hydrodynamic prediction for pr > 2GeV/ c |2l|. One 
might expect this to occur at lower px at ^/snn = 17GeV 
since the multiplicity is lower. 

NA44 is a focusing spectrometer EH- The dipole and 
quadrupole magnets produce a magnified image of the 
target in the tracking detectors. This has the advantage 
that particles that are close in momentum are not nec- 
essarily close in position. Since all tracking is done after 
the magnets we only have to reconstruct straight lines. 
A high resolution pad chamber sits in the focal plane and 
gives the magnitude of the momentum. Downstream of 
the focus the direction of tracks is measured by strip 
chambers and scintillating hodoscopes. The momentum 
resolution, including effects of multiple scattering in the 
target, is ~ HMeV/ c for p x , p y , and p z . 

A beam counter selected events for which single lead 
ions hit the 3.8mm lead target. Behind the target, two 
scintillator bars were used to trigger on high multiplic- 
ity events. An annular silicon detector with 512 pads 
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which was not part of the trigger measured dN^ jdr\ (r] = 
In tan 8/2) in the pseudorapidity range 1.5 < T) < 3.3. We 
used events in the top 18% of the multiplicity distribu- 
tion as we did for our earlier PbPb results 0, Q . For 
all data sets except the lower tut protons the spectrom- 
eter was set to accept particles of momentum 6-9 GeV/c 
and was positioned at 131 mrad with respect to the beam 
axis. The setting gave a pr window of 0.7-1.4 GeV/c. 
For the lower tut proton sample the momentum window 
was 5.2-8.0GeV/ c and the spectrometer angle was 44 
mrad. The resulting pt range was 0.18-0.50GeV/c. The 
rapidity range of the data was 2.4-2.9 and the system is 
symmetric about y=2.9 the center of mass rapidity. The 
Pt and y acceptance can be found in ■ 

Particles were identified by combining time of flight 
measurements from the hodoscopes with velocity infor- 
mation from three threshold-type gas Cherenkov coun- 
ters (CI, C2, TIC). C2 was set to fire only on pions. For 
the higher momentum proton data CI was used to re- 
ject kaons. For the lower m-r proton data kaons were re- 
jected exclusively by their flight time. For the tt~ sample 
the TIC was used to confirm that each track was a pion 
|23| . Contamination of the samples by other particles 
was less than 2%. However there is significant feed-down 
of hyperons into the proton sample, ranging from 28% 
for the lower m-r sample to 13% for (m-r) = 1-45 GeV 
|22|. This tends to smear the correlation function and 
was accounted for by including a fraction of "fake" pro- 
tons in the Monte Carlo 24] . The spectrum of hyperons 
was taken from RQMD with the yield scaled to match 
NA49 and WA97 results HHESl- 

The two-particle correlation function is defined by 



C 2 (pi,p 2 ) 



P 2 {pi,P 2 ) Real(p!,p 2 ) 



Pi(pi)Pi(p 2 ) Back(pi,p 2 y 



(1) 



where the numerator is the joint probability of detect- 
ing two particles with momenta p\ and p 2 , while the de- 
nominator is the product of the probabilities of detect- 
ing the single particles. The denominator was obtained 
by mixing tracks from two randomly selected different 
events. For protons, the attractive effect of the strong 
interaction competes with the negative correlation due 
to Fermi-Dirac statistics and Coulomb repulsion. This 
results in characteristic "dip-peak" structure in the two- 
proton correlation function. The height of the correlation 
peak decreases as the source size increases [27ll2^ . l29ll30j . 

The true correlation function is distorted by the mo- 
mentum resolution, Coulomb repulsion and by residual 
two particle correlations in the background. The magni- 
tude of these effects depends on the size of the source. 
For the pions we we evaluated them using an iterative 
Monte Carlo (MC) simulation. A trial correlation func- 
tion was used to generate events that we tracked through 
the spectrometer and analyzed in exactly the same way 
as the data. The Coulomb wave method was employed 
to simulate Coulomb effects |3l| . This method was used 



to derive correction factors to the correlation function 
and the process was repeated until the changes in the 
corrections from step to step were insignificant |lfl \?>2\. 
The theoretical correlation function for the pp data was 
generated by selecting proton pairs from a given source 
model, calculating the weight due to quantum statistics 
and final state interactions (strong and Coulomb) and 
propagating the particles through the Monte Caro. We 
searched for a source distribution that produced the min- 
imum y w ith the data after it was fed through the Monte 
Carlo :24j. 

C' 2 (pi,p 2 ) is a function of 6 variables; 3 for the total 
momentum P = p\ + p 2 , and 3 for the momentum dif- 
ference Q — p\ — p 2 . We decompose Q into Ql, the 
projection of the momentum difference onto the beam 
axis, Qout which is parallel to P and Q S ide which is per- 
pendicular to Qout and Ql- The data are analysed in the 
Longitudinal Co-Moving System in which p z \ — —p z2 . 

To measure correlation functions efficiently we trigger 
on events with pairs of particles that have a small mo- 
mentum difference in one direction while allowing large 
momentum differences in the other two dimensions. In 
our "vertical setting" the quadrupoles produce a small 
momentum acceptance in the horizontal space p x , and a 
wide acceptance in the vertical direction p y . Since the 
spectrometer lies in the horizontal plane this setting al- 
lows only a small range in Q S ide but a large range in 
Qout- In the "horizontal setting" we have a wide mo- 
mentum acceptance in p x and a small acceptance in p y . 
This setting allows only a small range in Q S ide but a 
large range in Q ou t- Both settings have a large accep- 
tance in Qhong- The lower ttit proton data were taken 
in the vertical setting while all other data were taken in 
the horizontal setting. 

For pions we used the maximum likelihood method to 
fit C2 with 



C 2 (Q T , Ql) = D(l + x e -Ql^-QlRl 



(2) 



where Rt and Rl parametrize the size of the source in 
the transverse and beam directions and D is a free pa- 
rameter used for normalization. The A factor is a mea- 
sure of the chaoticity of quantum states of the source; 
the fraction of pions from resonances and experimental 
factors that decrease the correlation function. For these 
data Qt ~ Qout- Figure ^ shows projections of the cor- 
relation function and the fit onto the Qt and Ql axis. 
The A parameter is larger than for our lower pt pion 
measurements. This is expected since the effect of reso- 
nances, which tend to wash out C 2 and reduce A, should 
be less important as pt increases. The systematic uncer- 
tainties in the fit parameters reflect the effect of (i) cut 
parameters to define a track, (ii) cut parameters to select 
pairs, (iii) momentum resolution, (iv) two-track resolu- 
tion, (v) momentum distribution of particle production 
in MC, and (vi) fitting to finite bins. 
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FIG. 1: The ir~ correlation function for (mj) — 0.92GeV pro- 
jected onto the (a) Qt and (b) the Qtaxies. Only statistical 
errors are shown. The lines are projections of the fit. 
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FIG. 2: Projections of the pp correlation function onto the 
Qinv axis for samples with different (mr). Only statistical 
errors are shown. The lines show the projections of the fit. 



Figure[3shows the pp correlation function as a function 
of Qinv = \/Q 2 — AE 2 for four bins in (tot) ■ At (tot) = 
1 GeV our result is consistent with NA49 I.331. Ri nv 
increases from pPb to SPb and finally to PbPb, |2J] . We 
have also analysed the pp correlation in two dimensions 
to extract Rt and Rl- Table [I] summarizes our data and 
compares them to two models. 

RQMD 2.04 is a string model which includes reinter- 
actions between produced hadrons (34|. It gives a good 
description of the proton data, as it does for our low tot 
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PbPb source radii versus hit from data and mod- 
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pion results 0, |3f| . The decrease of the radii with mj- 
has already been noticed in the model for pions and re- 
sults from correlations between position and momentum 
caused by flow j3ij. This decrease is more rapid for Rl 
than for Rt since the longitudinal flow is stronger. Nexus 
2.00 uses pertubative QCD to treat hard nucleon-nucleon 
interactions and pomeron exchange for soft interactions 
[3?|. It produces proton radii that are too small. Turn- 
ing on hadronic reinteractions and resonance production 
results in larger radii closer to the data. RQMD fails 
to reproduce both pion radii at tot = 0.92GeV. Nexus 
reproduces Rl but overestimates Rt by a factor of two. 

Under certain conditions [ljj one can derive analytic 
expressions for the radii: 



1 



1 



T^Af] 2 



1 ± A?7 : 



2%j 



(3) 



(4) 



where R is the geometrical size of the source, tq and To 
are the time and temperature of freeze-out, 8rj represents 
the width in rapidity and 77/ is the transverse rapidity. 
The radii also depend on the particle's mass but this 
effect is very small for reasonable values of 77/. Since 
r/f < At], Rl should drop more rapidly with ttit than 
Rt because the longitudinal flow is stronger than the 
transverse flow. 

Figure shows the m T dependence of Rt and Rl for 
all NA44 correlation and coalescence measurements 0, 
Hil l38| . The proton coalescence analysis gives only one 
radius which is mainly sensitive to the transverse size, 
(see Eqn. 6.3 of jl^|). These radii are in good agrement 
with the pp correlation data. Rl decreases faster with 
rriT than does Rt as expected from RQMD simulations 
and Eqns and The lines in Fig. [31 are the result of 
fitting our systematics to the form 1/R 2 = pO +pl ■ Mt- 
The bands show the errors on the fits. The x 2 /ndf is poor 
because the 7r~ radii at (tot) ~ 0.92 GeV are smaller 
than the trend expected from our other data. Removing 
the high itit tt" radii from the fit gives % 2 /ndf « 1 but 
has a negligible effect on the fit parameters. 
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FIG. 3: NA44 pion, kaon and proton source sizes versus hit 
[Hill HI- Both radii are fit to the form 1/R 2 = pO+pl-M T . 
The lines and bands show the results and errors of the fits. 



At tot = 0.22GeV the negative pions have source radii 
that are 9-11% smaller than for the positive pions. How- 
ever the 7T _ and n + correlation functions themselves are 
consistent and non Gaussian in shape 13]. RQMD simu- 
lations suggest that this is due to resonance decays, par- 
ticularly of the lu. Resonances tend to reduce the A pa- 
rameter when fit ting C*2 to Eqn. [21 This effect decreases 
rapidly with px |36| . It is possible that interaction with 
the residual nuclear charge could reduce the 7r~ radius 
as compared to the 7r + radius. Again this effect should 
be stronger at low pt- 

In summary, we have made the first two-dimensional 
measurement of the pp correlation function in relativis- 
tic nucleus-nucleus collisions. NA44 has measured HBT 
radii over a significantly larger tot range than any other 
experiment. The pp data, and all our previous HBT 
results, are consistent with a hydrodynamical interpre- 
tation and microscopic models that include hadronic 
rescattering and transverse expansion. Rl drops more 
rapidly with tot than Rt because the longitudinal flow 
is stronger than the transverse flow. At tot = 0.92GeV 
we find that the 7r _ radii are somewhat smaller than ex- 
pected from the trend of our other data and are not eas- 
ily explained by either hydrodynamic fits or the RQMD 
model. It may be that these fast pions left the source 
before the buildup of flow and so offer a glimpse of the 
hadronic system at an early point in its expansion. 
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